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PHASE-DIVERSITY WAVEFRONT SENSOR 

The present invention relates to a novel wavefiront sensor having a novel aberration 
means which can be constructed as a diffractive optical element (DOE) and which is 
suitable for use in adaptive optics. 

Adaptive optics (AO) is a technique in which real-time modification of the properties 
of an optical system based on measured optical output is used to indicate system 
optical performance. AO systems are widely used for compensation of turbulence in 
imaging through turbulent media in astronomical* and military** applications, for 
control of laser-beam delivery*", for control of laser resonator properties*^ for 
ophthalmic applications^ and in other areas. Li addition, aberration means can be used 
for measuring the distorting effect of a mirrored surface on an incident wavefiront 
when it is reflected firom the mirrored surface. 

AO systems typically comprise 3 components: 

• a wavefiront modulator (WFM) which alters the optical properties of the system in 
response to a command signal; 

• a wavefront sensor (WFS) which monitors the difference in state between the 
desired optical performance and the current optical performance of the system; 
and 

• a control loop which drives the WFM in response to the oulput firom the (WFS). 

For control of an adaptive optical system it is not necessary to reconstmct the input 
wavefiront. A sufficient condition for satisfactory operation of an adaptive optical 
system is the ability to drive a wavefiront modulator using a null sensor. In these 
circumstances, a control signal derived from a wavefront sensor system is used to 
indicate the size and, preferably, the location and the direction of the wavefront error. 
Thus, if the wavefront modulator is providing full correction of the input wavefront 
error, the control signal will be zero and the wavefront modulator will not be driven 
from its present position. The simplest example of such systems are based on multi- 
dither techniques applied, for example^*, to the correction of thermal lensing in laser 
resonators. 
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The phase-diversity^* algorithm can also be used to reconstruct the wavefront phase 
by using data corresponding to images of the input wavefront intensity on two planes 
normal to the direction of propagation and located at different positions along the axis 

5 of propagation. The approach taken is generally close to the 'two-defocus' method 
used in microscopy^" and operates with two close-to-focus images"^. Although the 
data planes are generally described as symmetrically placed about the image plane, 
they can equally well be symmetrically placed about the system input pupir, in which 
case the phase diversity algorithm becomes essentially the same as the wavefront 

10 curvature algorithm". 

For measurements on planes symmetrically spaced either side of either the image or 
the pupil planes, the intensity on the two measurement planes will be identical, and 
the difference between the images will be zero, if and only if, the wavefront in the 
15 . entrance pupil plane is a plane wavefront. 

The measurement of the difference between the intensity on the two data planes (the 
phase-div^e data), thus satisjBes the requirements for a null sensor because it 
provides a zero result when a plane wavefront occurs at the (input) pupil plane , If 
20 the input wavefront is distorted, the propagation between the measurement planes 
results in excess convergence (concave wavefront) or excess divergence (convex 
wavefront) and the resulting intensity difference between the measurement planes is 
indicative of the location, magnitude and direction of the wavefront curvature. 

25 Figure la shows a concave wavefront 3 propagating from left to right. For tiie 
sampled rays 3 shown, there is a local increase of intensity on one plane 5 and a 
decrease on title other plane 9. A convex wavefront 7 leads to a decrease in the local 
intensity. In each case, the intensities recorded in the +1 and -1 diffraction orders 
represent the intensities on the ends of the cylinder. 

30 

Figure lb shows that ambiguity can arise if the curvature of the wavefront 11 is so 
severe that a focus point occurs within the volume sampled by the measurements. 
The upper part of each figure leads to the same intensity increase. 
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Data on two image planes can be recorded using various approaches, including 
physical displacement of the image plane'', use of a vibrating spherically-distorted 
mirror, beam splitters and folded optical paths, or by the use of off-axis Fresnel 
lenses^^. Each of these approaches allows the recovery of phase information from two 
5 data sets recorded under different focus conditions. As such, all of these approaches 
are related to the two-defocus methods applied in electron microscopy in the 1970s. 

Altematively, a matched filter approach^", can be adopted by using a dif&active 
optical element to provide a pair of spots whose axial intensity (passing through a 
10 pinhole) defines the amplitude of each aberration mode in a distorted wavefront 

For reconstruction of the wavefront phase the inverse problem may be solved 
iteratively^'"^ or presented in terms of the differential Intensity Transport Equation 
(ITE) and solved through the use of Green's fimctions^^ with the phase-diverse data 
15 providing an estimate of the axial derivative of the intensity. Assumptions imposed 
through use of the ITE involve umformity of the input intensity in the entrance pupil, 
continuity in the wavefront phase and continuity of the first derivative of the 
wavefront phase'^^. 

20 At the expense of assumptions imposed on the uniformity of the input illumination, 
the work on phase-diversity has progressed a long way from the early electron 
microscope applications of the two-defocus technique and ITE-based approaches have 
demonstrated real-time data reduction with high (sub-nanometre) accuracy^^\ 
However, in all cases tihe two data sets are recorded under conditions where the 

25 wavefront is subject to a known defocus aberration between the two measurements. 

The restrictions on the uniformity and continuity of the input wavefront mean that 
current wavefront sensors cannot be used with scintillated, discontinuous, or multiply- 
connected wavefronts. There are many areas in which a generalised method which can 
30 cope with scintillated, discontinuous and disconnected wavefront segments would be 
a great advantage. One example is in polishing machines, where it would be useftil to 
measure in situ the surface profile of the object being polished. This would inevitably 
involve laser illumination of rough surfaces. It would also be usefiil to be able to 
image silicon circuitry, which by definition, is highly discontinuous. In addition. 
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modem large-aperture astrononaical telescopes are likely to be constructed using 
segmented, thus multiply-connected, mirror segments. 

It is an object of the present invention to provide a wavefront sensor which uses the 
5 null sensor concept that will produce a null output if the input wavefront is a plane 
wave, but which will retum an error signal when the input is distorted. A plane wave 
is defined as having constant phase across its wavefront. Without loss of generality 
that phase can be taken to be zero, so that its Fourier transform of the input 
wavefimction is Hermitian. 

10 

In accordance with a first aspect of tiie present invention there is provided a 
measuring apparatus for determining data relating to the shape of an input radiation 
wavefront, the wavefront shape being describable at a pre-determined location in an 
optical syst^: 

15 the apparatus comprising aberration means, the shape of which is defined by a filter 
fimction; 

detection means having a radiation sensitive surface capable of detecting the intensity 
of incident radiation on the surface, the detection means beiug coupled to an output 
device that provides a measure of the intensity of the incident radiation; 
20 wherein the aberration means acts on any input wavefront shape to produce first and 
second output radiation signals that in combination provide data from the output 
device on the extent to which the wavefront shape is non-planar. 



Preferably, the aberration means creates at least two filter Amotions, these filter 
25 fimction being complex conjugate pairs. 

Preferably, the filter fimctions are complex valued and have non-nuxed symmetry. 
This feature allows the invention to be used for the analysis of scintillated 
discontinuous or disconnected wavefront segments. 

30 

The filter fimction cannot be described by a pure defocus or an aberration means 
consisting of a diffractive optical element with essentially quadratic curvature in the . 
grating. Defocus diffractive optical elements cannot be used in the analysis of 
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scintillated, discontinuous or disconnected wavefront segments without generating 
boundary errors related to the wave shape. 

Preferably, the measuring apparatus further comprising a wavefront modulator. 

5 

Preferably, the wavefront modulator is configured to transform a wavefront being 
describable by means of a complex frinction into a wavefront being describable by a 
real ftmction. 

10 Preferably, the output device is provided with calculation means for calculating the 
difference between the first and second radiation signals. 

Preferably, the radiation sensitive surface of the detection means is provided with 
elements that allow the measurement of radiation intensity at different points across 
15 the siirface of the detection means. 

Preferably, the output device is provided with calculation means for calculating the 
difference between the first and second radiation signals at different points across the 
surface of the detection means. 

20 

Preferably, the wavefront modulator is coupled to the output device such that the 
wavefront modulator is distorted to provide a correction to a non-planar input 
radiation wavefront. 

25 Preferably, the aberration function is a weighted sum of Zemike polynomials arranged 
to equalise the signal generated from each mode of deformation in tiie input wavefront 
according to the expected statistical distribution of such modes in the input wavefront 

Preferably, the aberration means is arranged such that the complex conjugate filter 
30 functions of the aberration means are associated with diffraction orders of the same 
order but having different sign. 



Preferably, the first and second output radiation are produced simultaneously. 
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Optionally, the first and second output radiation are produced sequentially. 

Preferably, the aberration means is a diffractive optical element 

5 Optionally, the aberration means can be a variable-shape optical mirror 

Optionally, the aberration means is a variable refiractive index device. 

Optionally, the variable refiractive index device is a liquid crystal phase modulator 
10 used sequentially to provide complex conjugate aberrations or to provide a 
progranmiable diffiractive optical element. 

Optionally, the aberration means is a deformed reflective surface where the 
illumination of that surface firom each side produces the complex conjugate aberration 
15 fimctions. 

In accordance with a second aspect of the present invention there is provided a 
method for determining data relating to the shape of an input radiation wavefiront, the 
wavefi-ont shape being describable at a pre-detennined location in an optical system: 

20 the method comprising the steps of 

transmitting said input radiation wavefiront through an aberration means, the shape of 

which is defined by a filter fimction 

detecting fiie intensity of incident radiation on a surface, 

sending the detected intensity to an output device that provides a measure of the 

25 intensity of the incident radiation; 

wherein the aberration means acts on any input wavefiront shape to produce fibrst and 
second output radiation signals that in combination provide data from the output 
device on the extent to which the wavefront shape is non-planar. 

30 Preferably, the iterative reconstruction is used to reconstruct the wavefront shape. 



Preferably, the filter fimction is complex valued and has non-nuxed symmetry. 
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In accordance with a fhird aspect of the present invention th^e is provided a 
diflfractive optical element comprising abenration means, the shape of which is defined 
by a filter fimction that is complex valued and which has non-mixed symmetry. 

5 Preferably, the aberration means can be a variable-shape optical mirror 

Preferably, the aberration means is a variable refiractive index device. 

More preferably^ the variable refiractive index device is a liquid crystal phase 
10 modulator used sequentially to provide complex conjugate aberrations. 

Preferably, the aberration means is a deformed reflective surface where the 
illumination of that surface firom each side produces the complex conjugate aberration 
fimctions. 

15 

In accordance with a fourth aspect of the present invention there is provided a method 
for creating a dififractive optical element comprising the steps of defining aberration 
means the shape of which is defined by a filter fimction and which has non-mixed 
symmetry, applying said filter fimction to an optical element to create a diffractive 
20 optical element. 

Sufficient conditions for the null sensor are that it should give null output for input 
wavefimctions whose Fourier transforms have Hermitian symmetry, and that it give 
non-zero output for non-plane input wavefir^onts. In addition, the filter fimction should 
25 be complex. Hermitian symmetry exists where a fimction reflected in the origin is the 
complex conjugate of itself. 

Conditions where the real or imaginary part of the filter fimction is zero should be 
avoided to prevent a situation arising where the sensor will give null output for any 
30 input. In addition, mixed symmetries of the real and imaginary parts of the filter 

fimction should be avoided, since such filter fimctions provide an output signal even 
for a plane wave input 
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For imambiguous sensor output the real and imaginary parts of the filter function 
should have either both odd or both even symmetry. From these conditions it is 
possible to construct a null sensor without having placed any limiting assvimptions on 
the input wavefront. This method will therefore be suitable for scintillated and 
5 discontinuous input wavefironts. 

In accordance with a fifth aspect of the invention, there is provided a wavefiont 
sensor comprising: 
imaging means 
10 aberration means 
detecting means 
feedback means, 

arranged such that the imaging and detecting means capture images of the wave&ont 
field in a manner such tiiat the wavefiront field imaged is subject to two or more 
15 different aberrations imposed by the aberration means. 

Preferably, the aberration means when combined with the imaging means can deliver 
to the detecting means simultaneously or sequentially images of the input wavefi-ont 
field subject to different but known aberration point spread functions (psf) such that 
20 the filter functions (Fourier transforms of the aberration psfs) have complex odd or 
complex even symmetry and at least two of the images are subject to filter functions 
that are complex conjugate pairs. 

The aberration means does not provide pure defocus aberrations. 

25 

The filter functions may be described by any pure Zemike polynomial. 

Preferably, an electronic feedback means is arranged to deliver to an adaptive optics 
wavefiront modulator a signal derived firom the difference between the images affected 
30 by the complex conjugate filter functions and used to drive the wavefiront modulator 
in such a manner that the feedback signal strength is zero when the images recorded 
with complex conjugate filter functions are identical. 
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Preferably, the filter functions can be a weighted sum of Zernike polynomials 
arranged optimally to equalise the signal generated fi-om each mode of deformation to 
be corrected in the input wavefront field according to the expected statistical 
distribution of such modes in the input wavefi-ont field. 

5 

Preferably, the aberration means is a diffiractive optical element arranged such that the 
complex conjugate filter functions are associated with dif&action orders of the same 
order but different sign (±1, ±2, etc). Alternatively, the aberrating means can be a 
variable-shape optical mirror or a variable refi-active index device such as a liquid 
10 crystal phase modulator used sequentially to provide complex conjugate filter 

fimctions. Optionally, the aberration means can also consist of a deformed reflective 
surface where the illuxnination of that surface firom each side produces the complex 
conjugate filter fimctions. 

15 The following analysis show that the filter fimctions created by the aberration means 
are capable of encoding the sense and location of wavefront errors on heavily 
scintillated and discontinuous input wavefronts. The necessary and sufficient 
conditions required of filter fimctions suitable for use in a generalised phase-diversity 
null wavefront sensor for adaptive optics applications are quantified. 

20 

Let \P (r ) = (r)| e"^"^ represmt the complex amplitude distribution in the entrance 
pupil of an optical system, r being the co-ordinate in the pupil plane. If ^(r) 
represents a plane wavefront, the phase satisfies q> (r ) = constant . Any spatial 

variation in the wavefront phase represents a distortion, or aberration, that requires 
25 correction in an AO system. 

To operate successfully as a null wavefront sensor in an AO system we require a 
device that produces an error signal where, ^(r) is not constant. For the purposes of 

this analysis, the optical output from monochromatic systems are considered. Thus a 
30 discontinuity of exacfly an integral number of wavelengths in size does not affect the 
optical output, e.g. image intensity profile, and may be regarded as an undistorted 
wavefront. Such a discontinuity may be an important wavefront error in metrology 
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applications. It is desirable that the error signal provides information that localises the 
wavefront error in r-space and indicates the sense in which correction should be 
effected (in the absence of such indication a multi-dither technique^^" is reqxiired to 
effect correction). 

Let y^(^) = iy(^)+^(|) bethe Fourier transform of ^(r), where H{^) and A{^) 
represent respectively the Fourier transforms of the real and imaginary parts of ^(r). 
Clearly, H(^) is Hennitian and A{4) is anti-Hennitian. These symmetry properties 
will be required later. Thus, 

Let F± (^) = R (^) ± // (^) be tiie filter function that is the Fourier transform of an 
aberration psf (r) , which represents a complex function with which ^ (r) is 
convolved when forming an image of the system entrance pitpil. The functions R (^) 
15 and / ((^) are real-valued functions and the± indicates the use of two filter functions, 
in which the Fourier phase of the filter is reversed. In this work we are particularly 
interested in the necessary and sufficient conditions that constrain R (^) and /(^) in 

such a way that (^) are suitable filter functions to provide a null wavefront sensor 
for use in adaptive optics. 

20 

The physical relationship between these different functions is indicated schematically 
in figure 2 where wavefiront 1 1 is shown along with a diflBractive optical element 
which provides the aberration means 13. Images 15 and 17 with aberration function j 
+(r) and functions j-(r) respectively are shown along with an image 19 of the 
25 wavefiront. 



Hie detected intensity function may thus be written 



(2) 
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Substituting for , expanding and simplifying, the difference between the images 
fomied using the two filter functions may be expressed 
dir)^j\{r)-J.{r) 

= 2il^d^ V^(^)/(^)e-''« yr-C^O^COe*^*' (3) 

5 This is a real-valued function, since the quantity in [ ] is a difference of two complex 
conjugates and is thus imaginary-valued, so 

2i ^ ^ ' (4) 

The rhs of (4) can then be expanded and the terms grouped into 4 separate expressions 
10 which are equal to the rhs of (4) when summed: 

jd4Hi^)me-^^'jd^'H\4^^^^^ (5.1) 

Jj^ff(^)/(^K^^Jrf^'^\^')i?(^0^'"^ (5-2) 

15 jd^Ai^)Ii4)e-^^^jd4'H\4'm4'^^^^^^^ (5-3) 

jd^Ai^me'^^ldrA\^')Rine'^^''^ (5.4) 

This expression for the difference between the two detected intensity fbnctions is 
20 generally valid - no restricting assimiptions have so far been made. We may now 
investigate under what conditions of symmetry these expressions individually, or 
summed, are identically zcto. 



25 



Conditions for operation as a null wavefront sensor 
3.2.1 Filter function must be complex 
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Unless both R and / are non-zero, all of the terms in equations (5) are identically 
zero V I// . Thus d (r) is identically zero for all input wavefironts and no error signal is 
generated frona a non-flat wavefiront. 

5 This may be readily understood since if /(^ is zero is identical whatever the 
aritihmetic sign, so the two images are identical Vi/r . IfRiQ is zero the change of 
arithmetic sign is lost by the modulus square operation that is an inevitable part of 
quadrature detection processes — and thus the two images are again identical V \f/ . 

10 3 .2.2 Complex filter function 
3.2.2.1 Even symmetry 

Suppose that both /(^) and are even functions of ^ . 

15 

Consider the first expression (5.1). Since H[^) is Hemiitian and /(^) is symmetric 
and real-valued the product H{^)I(^ is Hermitian. Thus the Fourier integral 
ld^H{i^)I{^)e'^'^ is real-valued. The same is true of\d^'H\^')R{^'y'^^' . The 

second product of two integrals is term by term the complex conjugate of the first 
20 product. Thus expression (5.1), the difference between two complex conjugates, is 
always zero when bolli /(^) and R (^) are symmetric. 

Similarly (5.4) is always zero because each of the integrals reduces to a purely 
imaginary function. The product of these imaginary functions is real and the 
25 difference between the two complex conjugate terais is again always zero. 

Thus the difference between the two images is the sum of (5.2) and (5.3) and can be 
written 
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(6) 

In equation (6), for each pair of integrals one integral reduces to a real-valued function 
5 and the other integral to an imaginary-valued function. The first and second terms, 
also the third and fourth terms, are complex conjugate pairs and thus the rhs of (6) is 
imaginary valued or zero. 

If either H ot A is zero or if -4 = He^^ with q? constant, equation (6) and thus d(r) , 

10 is zero. However, these conditions are exactly those under which the input wavefiront 
is flat and a null sensor is required to produce a null output. 

Thus, if the filter function is complex with even symmetry the difference between the 
two images formed using these filters forms a potentially-useful null wavefiront 
15 sensor. 

3.2.2.2 Odd symmetry 

Suppose that both /(^) and R{^) are odd functions of ^ . 

20 

In (5.1) and in (5.4) the odd symmetry of the real functions R and / means that the 
products within the integrals in (5.1) have anti-Hermitian synmietry and those within 
(5.4) have Hermitian symmetry. Thus, although the arguments given are reversed for 
each term firom those given in 3.2.1, both (5.1) and (5.4) are identically zero Vy^ . 

25 Thus d (r) again reduces to equation (6). The arguments firom 3.2.1 again hold, 

although the role of the terms is reversed, one term in each integral product is purely 
real and the other term purely imaginary. 

Thus a filter function with complex odd symmetry is potentially suitable for use as a 
30 filter function for wavefiront sensing. 
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Suppose that one of /(^) and R (^) is an even function of ^ and the other is an odd 
functions of ^ . 

5 

In (5.1) the mixed symmetry will result in one of the integrals in each product being 
purely imaginary and the other integral purely real (dependent on whether / or i? is 
odd). In either case the product of the integrals is pvirely imaginary and thus (5.1) is 
purely imaginary or zero. An equivalent argument shows that (5.4) is purely 
10 imaginary or zero. 

Note, however, that the expressions (5.1) and (5.4) are reliant on the interactions 
between H and the filter function or on A and the filter function and do not involve 
cross terms between H and A . 

15 

In equation (6) the mixture of odd and even symmetry will mean that both integrals in 
each product are either purely real or purely imaginary. In each case the product of 
these terms will be purely real and the difference of the complex conjugates will thus 
be zero V^^. 

20 

Thus, filters with mixed symmetry produce no error signal dependent on the deviation 
of the wavefront from a plane wave and thus such filters are unsuitable for use as a 
wavefront sensor. 

25 3 .2,3 Sensing the error direction 

Information about the wavefront aberrations is encoded by the sum of the cross terms 
in equation (6). 

30 If the sense of the wavefront error is reversed, the phase of the wavefront will change 
arithmetic sign. This means that for a wavefront error of a given amplitude the error 
signal (equation (6) for filter functions with complex even or complex odd symmetry) 
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changes siga if tiie sense of the error is reversed. The wavefront sensor thus delivers 
an error signal that preserves information about the sense of the wavefront error. 

Note that, since the relationship between the error signal and the wavefront error is 
5 non-linear (involving the balance between H and A in the description of the 
wavefront), this does not guarantee that the error signal can be inverted to find the 
wavefront shape. 

3 ,2.4 Localising the wavefront error 

10 

Returning to equation (6) we note that each of the integrris, when expressed as a 
ftinction of r spears in the form of a convolution of ^ or -4 with a real-valued 
function deriving from the filter Amotion. As the products are of complex fimctions 
and real-valued ftmction the phase of the product is identically equal to that of the 
15 Hennitian or Anti-Hennitian ftinction in each integral. 

Without loss of generality, the location of the wavefront error can be identified with 
the position at which a(r) , the Fourier transform of A{^) , is non-zero. Because no 
phase ramp is added by the product with fhe real-valued filter fimction, the location of 
20 the non-zero componmt is, in some sense, localised at the point where a(r) is non- 
zero. 

The key term that encodes information about the wavefront aberrations is the sum of 
the cross terms in equation (6). 

25 

If the sense of the wavefront error is reversed the arithmetic sign of A[^) reverses 
because the phase of the wavefront will change arithmetic sign. This means tihat, for a 
wavefront error of given amplitude, the error signal (equation (6) for filter functions 
with complex even or complex odd symmetry) changes sign if the sense of the error is 
30 reversed. The wavefront sensor thus delivers an error signal that preserves 
information about the sense of the wavefront sensor. 
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Equation (6) also provides localisation of the wavefront error. Each of the integrals, 
when expressed as a function of r, appears in the form of a convolution of the Fourier 
transfonn of Hot A with a Hennitian function related to the filter function. Without 
loss of generality, the location of the wavefront error can be identified with the 
5 position at which a(r)y the Fourier transfonn of A (^) , is non-zero (a(r) is directly 

related to the phase distortion). The function witih which the convolution takes place 
may or may not have a maximum at the origin. When it does not the location of the 
non-zero component is, in some sense, localised at the point where a(r) is non-zero, in 
the sense that it may be distributed around the position at which a(r) is non-zero. 
10 When the intensities are recorded in pupil space this provides the localization required 
for a null wavefront sensor. 



The present invention will now be described by way of example only with reference 
to the accompanying drawings of which: 

15 

Figures la and lb show a series of wavefronts propagating and illustrates variations in 
local beam intensity; 

Figure 2 illustrates a wavefront and a diffractive optical element which provides the 
aberration means; 
20 Figure 3 illustrates an embodiment of the present invention; 

Figures 4a to 4d show simulation result with the non-plane wavefront input and the 
even symmetry phase diversity filter; 

Figure 5 shows a second embodiment of the present invention; 
Figure 6 shows a third embodiment of the present invention; 
25 Figure 7 shows a fourth embodiment of the present invention; 
Figure 8 shows a fiftii embodiment of the present invention; 

Figure 9 shows results with the Zemike polynomial Zio^ phase diversity filter used to 
create aberration means; 

Figures 10(a) (i) to (iii) show the wavefront intensity of a wavefront in the pupil 
30 phase, the output intensity distribution of a difference image and the associated prior 
art quadratic diffiractive optical element; 

Figures 1 1(b) (i) and (iii) show an embodiment of the present invention and illustrate 
the wavefix)nt intensity of a wavefront in the pupil phase, the output intensity 
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distribution of a difference image and the associated pure spherical diffractive optical 
element; and 

Figures 1(a) to 1 1(b) show a discontinuous pupil phase wavefront and outputs with a 
prior art defocus grating and a grating according to the present invention. 

5 

The requisite filter functions could be implemented in several ways, but one of the 
most straightforward approaches (for use with quasi-monochromatic radiation) is to 
use a diffractive optical element (DOE) to introduce the aberration. 

10 Figure 2 shows a DOE used to implement the requisite filter functions in the + 
dijBfraction orders and the zero order contains an image of the field to be 
reconstructed. In situations where the photometric efficiency is important a phase 
grating may be used and the depth of tiie grating adjusted to reduce to zero, or close to 
zero, the flux in tiie zero diffraction order. In Figure 2, the DOE 13 is shown that is 

15 used to generate the aberrations, the image with aberration psf j+(r), 1 5 and the image 
with aberration psf j.(r), 17 are also shown along with a representation of the 
wavefront to be reconstructed 19. From 19 we can tell whether the wavefront to be 
reconstructed is scintillated, multiply connected, uniform intensity etc. 

20 Figure 3 shows the combination of DOE 27 (diffraction grating) with SLM 29 (spatial 
line modiilator) operated using liquid crystals in this example. The wavefront 
modxilators are used to produce an adaptive optical system for driving a null sensor. 
In Figure 3 the CMOS (complementary metal oxide semiconductor) camera 23 used 
to detect the j+1 and j-1 signals is shown along with the lens 27and logic circuit 3 1 

25 which provides an electrical control system that is used to control the SLM 29 thweby 
altering flie shape of the wavefront in response to errors detected in the j+1 and j-1 
signals. 

In this example, the difference between the two images on a pixel by pixel basis is 
30 used directly to drive the wavefront modulators in order to provide the requisite 
wavefront correction. 



Spatial Light Modulators (SLM's) can be used to provide the wavefront modulation. 
SLM's provide a low-cost, low-voltage and low-volume modulator. In addition. 
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diversity can be applied to the input wavefront in the grating and lens combination. 
The grating is provided with a carefiilly selected filter function that satisfies the 
necessary and sufficient conditions for a null wavefiront sensor in accordance with the 
present invention. The grating allows the source image (in the 0th order) and the 
5 phase diverse data (in the ±1 orders) to be displayed on the focal plane, in which a 
CMOS camera is located. 

The images in the diffraction orders are subtracted to deduce the error signal and this 
information is fed back to the modulators to adjust the correction of the incoming 
10 wavefiront. 

A number of examples of the effect of a null wavefront sensor in accordance with the 
present invention are now provided. The above analysis has shown that a wavefront 
sensor can be created using a range of aberration means typically formed from 

15 dif&active optical elements. Furthermore, phase diversity filters can be used to give 
superior performance to the defocus method by using a priori information about the 
wavefront to be sensed. As seen in figure 2, image 19 shows the type of wavefront 
where the wavefront is scintillated, discontinuous or has discoimected wavefront 
segments. Aberration means is accordance with the present invention are used in the 

20 wavefront sensor to provide an accurate tool for determining wavefront error and 
therefore correcting the error in the wavefront. 

The embodiments of the present invention shown in Figures 3, and 5 to 8 of the 
present invention can be described as Compact adaptive optical systems (CAOS). In 

25 general, the apparatus is based on transparent wavefront modulators (e.g. liquid 

crystal SLMs) and phase diversity wavefront sensing. As described above. Compact 
Adaptive Optical (AO) systems, based on the use of DOEs, phase diverse wavefront 
sensing and transparent wavefront modulators, offer the potential for construction of 
AO systems with benefits such as: 

30 1 . Minimisation of non common-path errors by combining the wavefront 

sensor data and the corrected image in a single focal plane with essentially 
no separation of the 'science' and wavefront sensing optical trains; 
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2. The ability to place wavefront modulators conjugate to multiple planes in 
the object space whilst preserving a compact and robust optical train that 
would be particularly well-suited to multi-conjugate; and 

3. Avoidance of the requirement for conventional optics, such as beam 
5 splitters and conventional bulk optical components for re-imaging 

conjugate planes, that increases the size and weight of the equipment. 

Experimental simulations have been undertaken with plane and non-plane input 
wavefronts and a range of phase diversity filter functions having properties in 
10 accordance witii the present invention. 

Figures 4a to 4d show simulation result with the non-plane wavefiront input and the 
even symmetry phase diversity filter (Zio^, Zemike polynomial). From the simulation 
result we can see that the intensity differCTice between +1 order and -1 diffraction 
15 images is non-zero for the non-plane wavefiront input. This is a sufficient condition 
for a null wavefiront sensor derived from the theory. Figures 4 correspond to data 
recorded in an experimental embodiment of the form shown in figure 8. 

Figures 5 to 8 show additional embodiments of the present invention. 
20 For reasons of clarity, the detector/camera has been omitted from these figures. 

Figure 5 shows an embodiment 35 of the present invention in which an input/ test 
wavefit)nt 37 is shown at a position remote from a lens/grating 39. A modulator 41 is 
also shown connected to the output 51 of the logic drcuit 49 which is used to 
25 calculate the difiference in the intensity signals generated at Ii and I2 denoted by 

reference numerals 43 and 45 respectively. lo describes the position of the image of 
the wavefront 37. 

In this example of the present invention, the input/test wavefront is positioned 
30 remotely from the lens/grating 39. It has been foxmd that such an arrangement 

provides a sufficient degree of spatial resolution of the Ii and h signals to allow the 
actual position in the wave front of the non-planar part to be identified. This feature is 
particularly usefiil where the present invention is to be used as an adaptive optical 
device for the correction of input wavefronts . In this example, the measured , 
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spatially separated values of Ii and h will contain dark and bright areas where tihe 
wavefix)nt has deviated j&om being a planar wavefront. Once the difference between 
the calculated values of Ii and h has been calculated, the logic circuit 40 will create a 
signal which allows the modulator to adjust the shape of the wavefront to make it 
5 planar. 

Figure 6 shows a further embodiment 53 of the present invention in which llie 
input/test wavefront 55 is arranged in the plane of the lens/grating 39. In this 
example, the output at Ii and h provide less well spatially resolved signal. Therefore, 
it is more difficult to obtain information on the exact position across the wavefront 
10 where the wavefront is non-planar. In this example the logic drcuit is used to provide 
a good indication of the size of the error, but less spatial information. A value for the 
overall wavefront error may be obtained relatively quickly. Figures 4 and 9 show 
simulations of the data obtained for embodiments of the ^e shown in figure 6. 

15 Figures 7 and 8 provide additional embodim«its of the present invention with figure 7 
showing a wavefront remote from the lens/grating and figure 8 showing a wavefront 
in the plane of the lens/grating. In addition, a modulator 79 is provided in the 
apparatus of figure 8 to allow resolution of complex wavefront shapes. 

20 The examples of figures 7 and 8 are primarily intended for use with analytical 

instruments where the wavefront shape 59, 75 is a reflection firam a surfece sudi as a 
mirror. The reflections are used to measure the sh«tpe and quality of the surfece from 
which the wavefront has been reflected. As with the arrangement shown in figure 5, 
figure 7 provides an output from Ii and h that is relatively hi^y spatially resolved 

25 and which provides an output 51 in the form of data on the position of the wavefront 
error, which in turn provides information on the defect in the surface. As with the 
arrangement in figure 6, tiie arrangement in figure 8 provides a relatively quick 
measure of the magnitude of an error. 

30 Typically, these instruments can be used for measuring defects in astronomical 
mirrors. 
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Figure 9 shows results with the Zemike polynomial Zio® phase diversity filter 
function used to create aberration means in an embodiment of the present invention. 
This filter function is an even phase diversity. 

5 Figures 9a and 9b show the original signal phase and the retrieved signal phase. 
Figure 9c shows the difference between these phases, figure 9c shows the hi error 
value plotted against the number of iterations and figures 9e and 9f show the I+l and 
I-l images. 

10 Figures 1 0(a) (i) to (iii) show the wavefiront intensity of a wavefiront in the pupil 

phase, the output intensity distribution of a differmce image and the associated prior 
art quadratic diffractive optical element. Figures 10(b) (i) to (iii) show the wavefiront 
int^Qsity of a wavefiront in the pupil phase, the output intensity distribution of a 
difference image and the associated pure spherical dififractive optical element, an 

15 embodiment of the present invention. The prior art defocus grating is constructed 
using a curve described by the expression 2^ - 1. The spherical grating is 
constructed using a curve described by the expression 6^— €^ + 1 . 

The mathematical shape of the input wavefiront has sufficient continuity and 
20 tmiformity of intensity to operate with the prior art defocus grating and a grating 

suitable for use in the present invention.. As shown in the diagram of the pupil phase 
in Figures 10(a)(i) and 10(b)(i) are substantially identical. However, it will be noted 
that the amplitude of the input wavefront error for figure 10(b)(i) is approximately 
10% of the amplitude error shown in figure 10(a)(i). The difference images 10(a)(ii) 
25 1 0(b)(ii) have substantially, the same intensity. Therefore, the wavefiront sensor 

created using filter functions which are complex valued and which have non-mixed 
symmetry has superior sensitivity compared to the defocus only quadratic filter. 

Other non-quadratic polynomials may be used to define curves for constructing 
30 gratings in accordance with the present invention. 

The error graph is created using a known iteration algorithm and indicates that the 
error reduces to a substantially constant value after around SO iterations in this 
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example. It has been found that the error values for the Z\q polynomial are in this 
example, lower that those achieved using a defocus filter. 

Figures 11a and 1 lb show wavefront intensity in the pupil phase and output intensity 
5 after the wavefront has passed through a defocus gravity as shown in figure 10(a)(iii). 
Figures 11c and lid show wavefront intensity in the pupil phase and output intensity 
after the wavefront has passed througji a spherical aberration diffraction gravity as 
shown in figure 10(b)(iii). In both figures 1 1(a) and 1 1(c), the input wavefronts are 
identical and have a central obscuration 50 in the centre of the wavefix>nt. The central 
10 obscuration shown is of the type found in telescopes. 

It can be seen from figure 11(b) that there are a number of features around the 
boundary of the central obscuration 52 when the wavefront passes through a defocus 
gravity. These features arise because of the sensitivity of the defocus gravity to the 
15 imssing boundary and do not represent errors in the >yavefront, but are created by the 
optical properties of the defocus gravity. In a wavefront sensor, the presence of such 
features would prevent accurate modification of the wavefront shape and prevent the 
creation of a plane wavefront. 

20 It can be seen from figure 1 1 (d) that the boundary around the central obscuration 54 is 
smooth. Therefore, the use of a grating in accordance with the present invention is 
which the grating is constructed such that the shape of the aberrating means is defined 
by a filter fimction that is complex valued and has non-mixed symmetry. 

25 In fiirfher examples of the present invention, the input wave can be moved and the a 
grating whose shape is defined by a polynomial in accordance with the present 
invention, selected to provide optimum results. This can be achieved through 
statistical analysis of the wavefront. 

30 Improvements and modifications may be incorporated herein without deviating from 
the scope of the invention. 
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